Non-stationary magnetic microstructures in stellar thin accretion discs 
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We examine the morphology of magnetic structures in thin plasma accretion discs, generalizing 
a stationary ideal MHD model to the time-dependent visco-resistive case. Our analysis deals with 
small scale perturbations to a central dipole-like magnetic field, giving rise to the periodic modulation 
of magnetic flux surfaces along the radial direction, which corresponds to the formation of a sequence 
of toroidal current channels. These microstructures exhibit an exponential damping in time because 
of the non-zero resistivity coefficient, allowing us to define a lifetime of the configuration which 
mainly depends on the midplane temperature and on the length scale of the structure itself. If 
we try to identify a suitable phenomenological scenario where this framework could be applied, it 
turns out that the steady-state assumption adopted in literature is inconsistent with the evolution 
time-scale of typical stellar accretion disc systems; we can instead aim to describe peculiar transient 
events whose duration ranges from months or years. 
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I. INTRODUCTION 

A new perspective on the equilibrium morphology of 
thin accretion discs has been introduced in [TJ [2] , where 
it is outlined how the effect of the plasma back-reaction 
on the magnetic field of the central object can induce 
a radially oscillating crystal profile in the magnetic sur- 
faces (see also [3 0]), eventually breaking up the disc 
into a ring series, in the limit of a strong non-linear back- 
reaction. This issue is by its starting point different from 
the established picture of disc dynamics and its well- 
known open questions, like the problem of turbulence- 
enhanced viscosity (for review purposes, you can see the 
works of [5] and [5]). 

The brand new point of view stems from the imple- 
mentation of an ideal MHD scheme where the Ferraro 
Corotation Theorem [7j holds, allowing the angular ve- 
locity of the disc to be expressed only in terms of the 
magnetic flux surfaces (for a stationary visco-resistive ex- 
tension of this model see [SUH]). In fact, the background 
centrifugal force is balanced by the central object grav- 
ity - resulting in the Keplerian rotation of the disc - as 
well as by the background magnetic field, since the initial 
magnetosphere has a current-free morphology. When the 
back-reaction corrections to the centrifugal force and the 
Lorentz force are expressed in terms of the magnetic flux 
function, their balance results in the formation of toroidal 
current channels embedded in a radial oscillating mag- 
netic structure. A remarkable assumption in construct- 
ing this scenario is the pure rotation of the disc, which is 
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not endowed with any poloidal component of the veloc- 
ity field. Indeed, in many real astrophysical systems, the 
disc plasma is not very far from the quasi-ideal behaviour, 
but three main points are in order in this respect: 

(i) Given the mean thermodynamic parameters of the 
plasma in an accreting disc (temperature and number 
density), the kinetic theory uniquely determines finite 
non-zero values of the viscosity and resistivity coeffi- 
cients. Even when the value of such coefficients is rather 
small, the impact of their existence - i.e. the damping 
they can induce on the crystal structure - could be not 
negligible, due to the long lifetime of the accreting sys- 
tems. 

(ii) The Shakura model for accretion within a thin disc 
configuration [101 111) relies on very large values of the 
viscosity coefficient, able to balance the angular momen- 
tum transport responsible for a significant non-zero ac- 
cretion rate. Despite the crystal profile and the Shakura 
model of accretion are scenarios not directly comparable 
(e.g. in the crystal profile the purely rotating disc is un- 
able to directly accrete mass), it is of significant interest 
to understand if the periodic radial profile of the disc 
still survives in the presence of the visco-resistive effects 
required to account for the behaviour of a turbulent ac- 
creting plasma. 

(iii) A non-zero resistivity coefficient is always present in 
every numerical simulation, at least to reproduce the ef- 
fect of numerical dissipation due to machine's finite preci- 
sion [12] . None the less, works like [T3] include non-ideal 
terms to explore the coupling of disc and jet physics, 
aiming at the problem of angular momentum transport. 

In this work we include visco-resistive effects to the 
crystal profile, allowing for a time dependence of the mag- 
netic surfaces, needed to preserve the consistence of the 
Generalized Ohm Law. In fact, as far as we include a 
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resistive contribution, the azimuthal component of the 
Ohm Law acquires a term proportional to the toroidal 
current: this would turn out to be unbalanced in a purely 
rotating configuration, unless the non-stationarity of the 
model ensures the presence of a non-zero azimuthal elec- 
tric field. In such a non-stationary case, we are able 
to recover the Corotation Theorem and to provide the 
full consistence of the equilibrium configurations, using 
assumptions rather natural in the thin disc limit and re- 
quiring the toroidal component of the magnetic field to 
be negligible (we adopt a dipole configuration to charac- 
terize the magnetic properties of the central object). 

Focusing on the magnetic back-reaction, but neglecting 
its effect on the mass density distribution (note that this 
is the most general regime for a non-stationary purely 
rotating disc), we obtain a crystal profile for the disc 
magnetic field, isomorphic to the one studied in £Q, but 
for the presence of an exponential damping in time. The 
characteristic lifetime of this crystal structure depends 
on the resistivity (or equivalently on the viscosity, since 
the Prandtl number is constrained to be unity by the con- 
sistence of the configuration equations) and crucially on 
the spatial scale of the radial oscillations. Clearly, when 
such a lifetime parameter is greater than the average life 
of the considered accreting system, the resulting profile 
is mainly stationary and the analysis in pQ is recovered. 
However, a detailed study of the range of values of the 
damping lifetime leads to conclude that the crystal profile 
can not concern the stationary state of a stellar system, 
because the magnetic microstructures vanish on a time- 
scale much less then the typical lifetime of the accretion 
disc system. 

None the less, the obtained time-dependent scenario 
has the proper time-scale to account for cataclysmic phe- 
nomena. Their duration can be easily fitted in the range 
of density and temperature which are typical for stellar 
accretion discs, and in correspondence to radial scales 
that preserve the request to deal with a local model 
around a fixed value of the radius (i.e. the length scale 
of the oscillations is much below the disc radial size). 
The present analysis is therefore crucial in focusing the 
correct phenomenological scenario to which the crystal 
paradigm must be referred to in a stellar system, ruling 
out this morphology of the magnetic field from the steady 
state of an accreting stellar configuration. 



II. RELEVANT EQUATIONS 

Assuming electron pressure gradient is still negligible, 
the Generalized Ohm's Law retains its validity in the 
same form of stationary MHD: 



E+ ~[v- B ) = ;;./. 



(1) 



obviously meaningful at every different instant t. An 
important difference stands instead for the expression of 
electric field descending from Faraday's Law: 



E = -V$ - -d t A 



(2) 



written in terms of electrostatic potential $ and vectorial 
magnetic potential A. In particular, since the system is 
axially symmetric, we can say 
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where tp (r, z 2 ; i) is the magnetic flux function, leading to 

strictly poloidal components of the magnetic field B. 

Since we are dealing with a purely rotating configura- 
tion in local approximation (in a narrow annulus around 
a fixed radius), the azimuthal component of Generalized 
Ohm's Law (fil becomes: 



dtp rjc 2 f d 2 tp d 2 ip\ 
IH^Jtt [W 2 + W 2 = 



(4) 



Such an expression must be matched with the Induction 
Equation, namely 



— = -Vx [vxB - — -VxB 
dt \ 4tt 



(5) 



and this can be done only if the resistivity it's nearly a 
constant, precisely if 
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where A(-) is the Laplace operator. This condition ensure 
that radial and vertical components of Equation ^ give 
again Equation Q, so they are redundant, whereas the 
azimuthal one states the following constraint: 



VwxVV> = 0, 



(7) 



telling us that the angular velocity is a flux function 
Lj(r,z]t) — ui{ip). This result generalizes the Corotation 
Theorem of [7], still holding in the visco-resistive frame- 
work provided there is no azimuthal component of the 
central magnetic field. 

Finally, we deal with the MHD momentum conserva- 
tion equation: 



p(d t v + (v-Vj vj = - Vp - pV\ + Fl- 



+ v 



w 2 v + ^v(w 



(8) 



where p is thermodynamic pressure, x is t ne gravitational 
potential of the central object, Fl is the Lorentz force, 
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and the viscosity v is assumed to be a constant. Its 
azimuthal component gives us the evolution law for the 
angular velocity: 



~dt 



P 



d 2 io d 2 uj 
dr 2 dz 2 



0. 



(9) 



while the other components retain their stationary forms 
because of the assumption of pure rotation. 



III. PERTURB ATIVE SCHEME 

Now we can separate the contribution of the back- 
ground dipole-like magnetic field from the back-reaction 
induced by the plasma current, i.e. 



ip(r,z 2 ;t) =tPo(R ,z 2 ) + ^(r - R , z 2 ;t) 



(10) 



where Rq is the fiducial value for our local study, cen- 
tred on | r — Rq I <C Rq ■ We highlight that background 
■00 is a stationary vacuum solution of Laplace's equation 
(i.e. AtpQ = 0), so Equation Q can be written by means 
of back-reaction ipi only. 

Correspondingly, since Equation ^ constrains the an- 
gular velocity to depend upon the magnetic surfaces (and 
mainly on the back-ground one), an expression for the lo- 
cal study around Rq is the following: 



w(i? , r - Rq, z; t) = u(<ip) ~ u>(ipo) + oj' ipi , (11) 



with 
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Note now that the evolution of angular velocity is de- 
termined by Equation Q: using the expansion (111, it 
reduces itself exactly to Equation Q if 
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which states that the Magnetic Prandtl Number is of or- 
der unity, condition also adopted in [5] with the purpose 
of recovering the existence of a solution consistent with 
the Corotation Theorem. 



Furthermore, Equation ( 13 1 tells us that the height- 
dependence of r\ = rj{p) is determined only by the density 
profile, giving a deeper meaning to Equation This 
can be explained after the introduction of dimensionless 
variables defined as follows: 



x = k(r — Rq) , u 



^zHo 
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T PO 



and dimensionless functions written as 

2 A\ 
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Po 
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where po = p(x = 0, u = 0) is the density value on the 
equatorial plane at the fixed radius, Hq is the half- 
thickness of the disc, and A: -1 is the radial scale of the 
back- reaction flux function ip\. The parameter e z used 
in Equation ( 14 ) to scale the height-dependence of mag- 



netic surfaces is related to the plasma /3-parameter by 
the relation 
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with po the back-ground thermodynamic pressure. We 
can then write down the components of Equation Q to 
obtain 



drV 1 d r ipi 

which are identically satisfied if 
e z < 1, 



since Equation (17b) can be restated as 
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(17b) 
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and this is assured for a thin disc. Then it follows that 



Equation ( 18 1 specifies the regime of validity for our 
treatment. 

In terms of these variables and parameters, Equation 
Q becomes 
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(20) 



which is a linear parabolic PDE with infinite many solu- 
tions. In particular, it admits a separable solution ana- 
logue to those found in [5], currently endowed with an 
explicit time-dependence: 



Y(x,u 2 ;6) = AF(u 2 ) sin(ax)e" 



(21) 



with A and a some real constants, whereas the height- 
dependence is fixed by 



d 2 F (u 2 ) 

e z ^l-(a 2 -D{u 2 ))F(u 2 ) 



0, 



(22) 



meaning that it is ultimately fixed after the choice of an 
equation of state, needed to set an expression for D(u 2 ). 
If we choose a polytropic form for the gravothermal back- 
ground, such that: 



piu 



= PqD(u 2 ) 
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FIG. 1. Behaviour of the back-reaction magnetic surfaces 
in Equation (24 1 against dimensionless vertical coordinate, 
normalized to the gaussian equatorial value obtained at u = 
with a = a m ; n . We fix e z = 0.05 and T = 1, then we compute 
a in the range [a m ; n , 1.2]. As a increases, the plots are drawn 
as follows: thick, thin, dashed, dotted, dot-dashed. It is worth 
noting that the gaussian profile has the greater FWHM, while 
it hasn't the greater initial amplitude. 



where T is the polytropic index, then we can find an an- 
alytic solution of Equation ( 22 ) by means of Generalized 
Hcrmite Polynomials: 



F(u 2 ) 



1 



2e z Vf 



'Tu 2 



(24) 

having omitted a constant factor and adopted the no- 
tation H[£, x] for the Hermite Polynomial of order £ in 
the variable x. Note that these Generalized Polyno- 
mials are monotonically decreasing provided their or- 
der is negative. Choosing the radial wave-number as 
a 2 = a 2 lin = 1 — e z VT, the solution p4| reproduces the 
gaussian profile already obtained in |1J. This is the lower 
limit for the range of values of the parameter, since for 
a < flmin the solution becomes non-physical; an upper 
limit doesn't exist instead, although the amplitude of the 
structure becomes negligible as a increases (see Fig. (fl])). 
Therefore, only a narrow range of values around the unity 
is suitable to tune the radial wavenumber. 



tells us that the density has to keep its steady-state pro- 
file, which assures function D(u 2 ) doesn't depend on 
time. 

Radial and vertical components of momentum conser- 
vation ([8 ) give their stationary formulation too, returning 
exactly the same dimensionless equations offered in 
now coupled to Equation ( 20 1 . This equation introduces 
the new variable 9, but acts as a closure condition for the 
system which determines the equilibrium of perturbative 
pressure P and density D, namely: 

d u 2P + e z D + 2 (d 2 x Y + e z 8 2 u Y) d^Y = (26a) 



-d x P 



3eiD) Y+ 



+ (d 2 .Y + e z dlY)(l + d x Y) = 0, 



(26b) 



In the linear regime Y <C 1, when pressure and density 
perturbations are assumed to be negligible, system ( 26 ) 
reduces itself only to 



(dlY 



eMY) = -D( 



Y, 



(27) 



which is consistent with separable solution (21 1 and owns 
the same magnetic structure (although damped in time) 
developed in [5] or [5], according to the additional ap- 
proximations considered in those papers. It is interesting 
to note that Equation (20 1 could be a closure condition 



also in the general non-linear regime, but this is forbidden 
by Equation (25 1: this suggests that a complete model 



with the disc decomposition in a ring-like structure re- 
sembling the one in [2] (which needs the density to be 
time-dependent) should exhibit poloidal velocities. 

Another solution can be found in the special 'crystal 
regime' specified by 



Y>1, D<^D. 



(28) 



occurring when magnetic perturbation are relevant (i.e. 
back-reaction field is strong) whereas density perturba- 
tions are negligible. In this case the system ( 26 ) reduces 
itself to: 



d B Y - A F Y = 



d u iP + 2A £ Yd u 2Y = 



(29a) 



(29b) 



IV. MATCHING WITH PREVIOUS MODELS 



-d x P + D(u 2 ) Y + A e Y(l + d x Y) = , 



(29c) 



We underline that, under the assumption of pure rota- 
tion v = u;i?o60, the only equations explicitly depending 
on time are those shown in Section [TTJ although there are 
other interesting equations. Continuity equation in the 
form: 

dp 



dt 



+ V-(p0) = O, 



(25) 



which admits the solution ( 21 ) and gives a separable form 
to pressure too, explicitly: 



P(x,u 2 ;8) = A 2 D(u 2 ) F 2 (u 2 ) sin 2 (ax)e 



-20 



(30) 



where D[u 2 ) and F(u 2 ) can be expressed via Equation 
(23b) and Equation (24), respectively. 
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FIG. 2. Plot of microstructures' lifetime against their length 
scale via Equation (|3lJ for T e = 10 3 K (dashed), 10 5 K (dot- 
ted), 10 r K (dot-dashed), 10 9 K (thin line). Shaded region 
highlights lifetime values whose order of magnitude is from 
month to year. Thick horizontal line marks the value of 10 7 
years. See text for details. 



V. MEANING OF DAMPING TIME 

It turns out that every crystal magnetic structure pre- 
viously obtained as an equilibrium configuration is cur- 
rently a dynamical solution, after the compatibility with 
Equation (20) has been checked. At the same time, it 



has to suffer a damping like solution ( 21 ), which restores 



the unperturbed magnetic configuration dominated by 
the central field. Back-reaction microstructured field be- 
comes negligible after the time 



po 



(31) 



defined in Equation ( 14 ) . A direct consequence is that 



such structures could exist only if the plasma is quasi- 
ideal, with very low viscosity coefficient; in particular, 
this shows how the present analysis thwarts the Standard 
Model [ini E] ■ Let us try to use the phenomenological 
Shakura prescription for the viscosity, i.e. 



ap v s H , 



(32) 



where v$ is the sound velocity and a is a model parame- 
ter which determines the strength of turbulence enhance- 
ment. Remembering the Standard Disc is nearly Kep- 
lerian with fix — vs/Hq, and using Equation (32) in 



Equation (31), we obtain 



^standard 



-CI 



K ' 



(33) 



whose rough estimates (with the common choice of a 
in the range 0.01 — 0.1) are of the order of hours for 



a X-Ray Binary System, absolutely meaningless on the 
evolution time-scale. This is to prove how the scheme of 
magnetic crystal microstructures shows itself as an aut 
aut alternative to the Standard Model, since those kind of 
configurations are not held by the inbuilt assumptions of 
the a paradigm, at least if we seek for a quasi-stationary 
picture on the time-scale of stellar evolution. 



On the other hand, relation (31 ) provides us with con 



straints on the kind of systems in which those structures 
could be found. To this purpose, we use Equation (13) 
and the formula for microscopical resistivity: 



m e v ie 
n e e 2 



1.435 • 10" 



log A 



(T[K]) 



3/2 



(34) 



where the Coulomb Logarithm can be approximated by 
log A ~ 20 at the typical densities of interest. The mean 
lifetime r can then be expressed as a function of only 
the temperature T and the back-reaction length scale 
A = 271-fc- 1 . 

In Fig. ([2| we show r(A) for some different values of 
temperature. The horizontal thick line marks the esti- 
mated lifetime (10 7 years) for an accretion disc in High 
Mass X-Ray Binary Systems, when the disc can contin- 
uously replenish matter from the companion star. Mi- 
crostructures' lifetime has to exceed this value in order 
to speak of a quasi-stationary configuration on the scale 
of stellar evolution, but you can see that this threshold is 
achieved only with very high temperatures: a reasonable 
value of 10 7 K yields the length scale over 10 9 cm, which 
is comparable to the disc size; reducing this scale to a 
more meaningful 10 7 cm needs the temperature to ex- 
ceed 10 9 K, far beyond the typical range. We deduce this 
model is inconsistent with those time-scales and should 
not be applied in this way. We point out now that in [3] 
an estimation of A in the case of a global structure is pro- 
vided, and it is of the order of 10 3 — 10 4 cm, which gives 
rise to a negligible lifetime regardless of temperature. 



VI. CONCLUSIONS 

We discussed how the equilibrium configurations pro- 
posed by pQ and following works have to be generalized 
by means of non-ideal terms (an arbitrarily small but 
non-zero resistivity is required by kinetic theory, and 
its effects are not negligible), and how these equilibria 
evolve in time. The peculiar morphology of magnetic 
microstructures is preserved, but it has been shown how 
the dissipative coefficients forbids its maintaining: in par- 
ticular, the macroscopical viscosity of Shakura Standard 
Disc destroys it in a few hours or less, and the microscop- 
ical plasma resistivity is nevertheless unable to reach the 
time-scales of stellar evolution in binary systems. In spite 
of this, we noted that ranges of length and temperature 
exist (shaded region in Fig. ([2])) which are consistent 
with an evolution in terms of transient phenomena whose 
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lifetime is between months and years. Such a framework 
could involve plasma discs in the so called Cataclysmic 
Variables (e.g., DQ Herculis systems: see [2]), which are 
known to exhibit short-period events, as well as in Black 
Hole X-Ray Transient systems, which undergo very long 
periods of quiescence interrupted by non-periodic small- 
duration times of strong activity - for instance, see the 



atypical X-ray source GRO J1655-40 [15]. Temperatures 
in 10 3 — 10 5 K and length scales in 10 6 — 10 s cm are suit- 
able for this kind of objects, suggesting the formation and 
disruption of magnetic microstructures; the same range is 
available for the description of XRB systems in the case 
of no replenishing of matter, since isolated discs have 
similar lifetimes. 
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